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Abstract 

^^ ', The diffuse gamma-ray background radiation (GBR) at high Galactic latitudes could be 

o ■ 

(^ ' dominated by inverse Compton scattering (ICS) of cosmic ray (CR) electrons on the cosmic 

,_^ ■ microwave background radiation (GBR) and on starlight (SL) in an extended CR halo of 

^: 

'"^ , our galaxy. Assuming the locally observed GR electron spectrum beyond a few GeV, 

(N ■ 

^SJ I which follows from equilibrium between Fermi acceleration and radiative cooling, to be 

T-H ■ a universal spectrum throughout our galaxy and its CR halo, we reproduce the observed 

> ; 

0^ , spectral index, the intensity and the angular dependence of the GBR, in directions away 

"^ I from the galactic disk and centre, without recourse to hypothetical extragalactic sources. 

o: 

f—N I 1. The GBR: The existence of an isotropic, diffuse extragalactic gamma background 

i-Q ' radiation (GBR) was first suggested by data from the SAS 2 satellite (Thompson & Fichtel 

Q . 1982) and claimed to be confirmed by the EGRET instrument on the Compton Gamma 

H ; 

c/3 I Ray Observatory (e.g. Sreekumar et al. 1998). The alleged extragalactic diffuse GBR is the 

a . 

L^ ' diffuse emission observed by EGRET by masking the galactic plane at latitudes |b| < 10°, 

S^ . as well as the galactic centre at |bj < 30° for longitudes |1| < 40°, and by extrapolating to 

P^ \ zero column density, to eliminate the n^ and bremsstrahlung contributions to the observed 

radiation and to tame the model-dependence of the results. Outside the mask, the GBR 
flux integrated over all directions in the observed energy range of 30 MeV to 120 GeV, 
shown in Fig. 1, is well described by a power law: 



dF 

—^ ~ (2.74 ±0.11) X 10"^ 
dE 



E 



MeV 



-2.10±0.03 



cm~^ s"' sr"' MeV"' . (1) 



The overall magnitude in Eq. (|l]) is sensitive to the model used to subtract the foreground 
(Sreekumar et al. 1998; Strong et al. 1998), but the spectral index is not. The EGRET data 
are given in Sreekumar et al. (1998) for 36 (b, 1) domains, 9 values for each half-hemisphere. 



The spectral index is, within errors, extremely directionally uniform, as shown in Fig. 2, 
where we have plotted the EGRET results as functions of 9, the observation angle relative 
to the direction to the galactic centre (cos^ = cos [b] cos[l]). The normalization is less 
homogeneous, but in directions well above the galactic disk and away from the galactic- 
centre region it has been found to be consistent with a normal distribution around the 
mean value: thus the claim of a possible extragalactic origin (Sreekumar et al. 1998). 

The origin of this diffuse GBR is still unknown. The published candidate sources range 
from the quite conventional to the decisively speculative. Perhaps the most conservative 
hypothesis is that the GBR is extragalactic, and originates from active galaxies (Bignami 
et al. 1979; Kazanas & Protheroe 1983; Stecker & Salamon 1996). The fact that blazars 
have a 7-ray spectrum with an average index 2.15 it 0.04, compatible with that of the 
GBR, supports this hypothesis (Chiang & Mukerjee 1998). The possibility has also been 
discussed that Geminga-type pulsars, expelled into the galactic halo by asymmetric super- 
nova explosions, be abundant enough to explain the GBR (Dixon et al. 1998; Hartmann 
1995). More exotic hypotheses include a baryon-symmetric universe (Stecker et al. 1971), 
now excluded (Cohen et al. 1998), primordial black hole evaporation (Page & Hawking 
1976; Hawking 1977), supermassive black holes formed at very high redshift (Gnedin & 
Ostriker 1992), annihilation of weakly interactive big-bang remnants (Silk & Srednicki 
1984; Rudaz & Stecker 1991), and a long etc. 

The EGRET GBR data in directions above the galactic disk and centre show a signif- 
icant deviation from isotropy, correlated with the structure of our galaxy and our position 
relative to its centre (Dar et al. 1999). In Fig. 4 we have plotted, as a function of 9, 
the EGRET GBR counting-rate above 100 MeV. This figure clearly shows, in three out 
of the four quarters of the celestial sphere, an increase of the counting rate towards the 
galactic centre. The "reduced" x^ per degree of freedom for a constant flux is 2.6: very 
unsatisfactory. A best fit of the form F = Fq + Fi (1 — cos 9) yields x^ = 1.3, a very large 
amelioration (for higher polynomials in cos 9 the higher-order coefficients are compatible 
with zero: the fit does not significantly improve). Note also that at angles with cos 9 larger 



than its mean value {cos 6) = 0.0246 (6 < 88.6°), 10 out of the 12 data points are above 
the average flux, while at angles with 6 > 88.6°, 18 out of the 24 data points are below 
the average. The probability for a uniform distribution to produce this large or a larger 
fluctuation is 1.5 x 10^^. Even in directions pointing to the galactic disk and the galactic 
centre, EGRET data on 7-rays above 1 GeV show an excess over the expectation from 
galactic cosmic-ray production of vr^'s (Pohl & Esposito 1998). Electron bremsstrahlung 
in gas is not the source of the 1-30 MeV inner-Galaxy 7-rays observed by COMPTEL 
(Strong et al. 1997), since their galactic latitude distribution is broader than that of the 
gas. These flndings also imply that inverse Compton scattering may be much more impor- 
tant than previously believed (Strong & Moskalenko 1998; Moskalenko and Strong, 2000; 
Dar et al. 1999). 

Thus, it appears that the EGRET data advocates a local (as opposed to cosmological) 
origin for the GBR. Indications of a large galactic contribution to the GBR at large 
latitudes were independently found by Dixon et al. (1998) by means of a wavelet-based 
"non-parametric" approach that makes no reference to a particular model. Strong & 
Moskalenko (1998) and Moskalenko & Strong (2000) also found that the contribution of 
inverse Compton scattering of galactic cosmic ray electrons to the diffuse 7-ray background 
is presumably much larger than previously thought. 

2. Galactic Origin of the GBR: In recent publications we went one step further 
(Dar et al. 1999; Dar & De Rjula 2000) and explored in detail the possibility that the 
diffuse gamma-ray background radiation at high galactic latitudes could be dominated 
by inverse Compton scattering of cosmic ray (CR) electrons on the cosmic microwave 
background radiation and on starlight from our own galaxy. We made three very simple 
assumptions: (1) The locally-measured electron spectrum, shown in Fig. 4, which is well 
fitted, from E ~ 10 GeV to ~ 2 TeV by. 
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is representative of its average form throughout the Galaxy and its extended CR halo. (2) 
Above a few GeV the acceleration of CR electrons is in equilibrium with their cooling by 



inverse Conipton scattering (ICS) on starlight and on the microwave background radiation. 
(3) The GBR is produced by ICS of CR electrons on CBR and starlight photons in our 
Galaxy and in external galaxies. 

Assumptions (1) and (2) allow us to derive directly, the GBR spectral index from the 
observed electron index: (5^ = (/3c + l)/2 = 2.10 ± 0.05. It agrees with the GBR index as 
observed by EGRET, /3^ = 2.10 ± 0.03. This index is universal, independent of whether 
the ICS took place in our Galaxy or in external galaxies. 

The two dominant contributions to the GBR within our model are inverse Compton 
scattering of Galactic CR-electrons off the cosmic background radiation and starlight. 
There is a small additive effect of ICS from sunlight. The contribution from external 
galaxies, is also sub-dominant. The ICS spectrum -a cumbersome convolution (Felten 
&; Morrison 1966) of a CR power spectrum duc/dE = AE~'^° with a photon thermal 
distribution- can be approximated very simply (Dar & De Riijula 2000) by 
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where ctt = 0.65 x 10^^^ cm~^ is the Thomson cross section, ei are the mean photon 
energies of starlight and of the CBR and Ni(b, 1) are the effective CR electron column 
density resulting from the convolution of the space distribution of CR electrons with those 
of starlight and of the CBR in the direction (b,l). Simple analytical expressions for them 
are given in Dar & De Riijula (2000). 

Our predictions of the intensity and the angular dependence of the GBR are shown in 
Figs. 5, 6. Like that for the spectral index, they are in excellent agreement with those 
observed for the alleged extragalactic GBR, leading to the conclusion that the GBR could 
be dominated by the emission from our own galaxy. 

The only non-conventional aspect of our model is that, in order to reproduce the 
observed intensity of the GBR, we must assume the scale height of our galaxy's CR- 
electron distribution to be almost twice the traditionally-accepted upper limit. Such a 



large scale height is not in contradiction with radio synchrotron-emission from our galaxy 
if the galactic disk and its magnetic field are embedded in a larger magnetic halo with a 
much weaker field. 

3. Predictions: The main predictions specific to our scenario are: 

(1) The GBR should reflect the asymmetry of our off-centre position in the Galaxy. 

(2) Very nearby star-burst Galaxies, such as M82, and radio galaxies with presumably 
large CR production rates, such as Cygnus A, may be visible in gamma rays. 

(3) The GBR spectrum should not have the sharp cutoff, above E ~ 100 GeV, expected 
(Salamon & Stecker 1998) for cosmological sources. But it should nonetheless steepen 
around 10-100 GeV, because of the anticipated "knee" in the electron spectrum and of 
the energy-dependence of the Klein-Nishina cross section. 



These features of our scenario should be testable when the next generation of cosmic-ray 
and 7-ray satellites (AMS-02 and GLAST) are operational, hopefully by 2005. 
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Figure 1: Comparison between the spectrum of the GBR, measured by EGRET (Sreeku- 
mar et al. 1998), and the prediction for ICS of starlight and the CMB by CR electrons. 
The slope is our central prediction, the normalization is the one obtained for he = 20 kpc, 
Pe = 35 kpc. 
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Figure 2: EGRET data on the GBR spectral index as a function of 9 the angle away 
from the direction of the galactic center. The line is the predicted spectral index. The 
various plots correspond to the individual half- hemispheres, (a) b > 0, 1 > 0. (b) b > 0, 
1 < 0. (c) b < 0, 1 > 0. (d) b < 0, 1 < 0. 
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Figure 3: EGRET data, organized as in Fig. 2, for the dependence on 9 of the GBR 
intensity above 100 MeV. 
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Figure 4: The primary cosmic-ray electron spectrum (Evenson & Meyers 1984; Golden et 
al. 1994; Ferrando et al. 1996) as measured by Prince 1979 [crosses]; Nishimura et al. 1980 
[squares]; Tang 1984 [circles]; Golden et al. 1984 [triangles]; Barwick et al. 1998 [stars]. 
The slope is the prediction, the magnitude is normalized to the data. 
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Figure 5: Contributions to the GBR flux above 100 MeV as functions of longitude 1, at 
fixed latitude b, from ICS of starlight (dotted), and CMB in our galaxy (dashed); from 
the total ICS from external galaxies (continuous), and from sunlight (dot-dashed). The 
vertical scale is 10'^ times the number of photons/ [cm^ ssr]. The results are for he = 20 
kpc, pe = 35 kpc. 
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Figure 6: The flux of GBR pliotons above 100 MeV: comparison between EGRET data 
and our model for he = 20 kpc, pc = 35 kpc, as functions of latitude at various fixed 
longitudes. The grey domain is EGRET's mask. 
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Abstract 

The diffuse gamma-ray background radiation (GBR) at high Galactic latitudes could be 

dominated by inverse Compton scattering (ICS) of cosmic ray (CR) electrons on the cosmic 

O 

(^ ' microwave background radiation (GBR) and on starlight (SL) in an extended CR halo of 

(N ■ 

our galaxy. Assuming the locally observed GR electron spectrum beyond a few GeV, 

which follows from equilibrium between Fermi acceleration and radiative cooling, to be 

(N 

(NJ I a universal spectrum throughout our galaxy and its CR halo, we reproduce the observed 

spectral index, the intensity and the angular dependence of the GBR, in directions away 
> , 
0^ , from the galactic disk and centre, without recourse to hypothetical extragalactic sources. 

^-^ ■ 1. The GBR: The existence of an isotropic, diffuse extragalactic gamma background 

radiation (GBR) was first suggested by data from the SAS 2 satellite (Thompson & Fichtel 



i-Q ' 1982) and claimed to be confirmed by the EGRET instrument on the Compton Gamma 

Q . Ray Observatory (e.g. Sreekumar et al. 1998). The alleged extragalactic diffuse GBR is the 

H ; 

c/3 I diffuse emission observed by EGRET by masking the galactic plane at latitudes |b| < 10°, 

a . 

as well as the galactic centre at |b| < 30° for longitudes |1| < 40°, and by extrapolating to 
S^ . zero column density, to eliminate the vr'^ and bremsstrahlung contributions to the observed 

5^ , radiation and to tame the model-dependence of the results. Outside the mask, the GBR 

flux integrated over all directions in the observed energy range of 30 MeV to 120 GeV, 

shown in Fig. 1, is well described by a power law: 

E 



jp r -o 1 ~2.10±0.03 

■^ -- (2.74 ±0.11) X 10"^ 



dE 



MeV 



cm-2 s'^ sr-i MeV"^ . (11 



The overall magnitude in Eq. (|l]) is sensitive to the model used to subtract the foreground 
(Sreekumar et al. 1998; Strong et al. 1998), but the spectral index is not. The EGRET data 
are given in Sreekumar et al. (1998) for 36 (b, 1) domains, 9 values for each half- hemisphere. 
The spectral index is, within errors, extremely directionally uniform, as shown in Fig. 2, 



where we have plotted the EGRET results as functions of 9, the observation angle relative 
to the direction to the galactic centre (cos^ = cos [b] cos[l]). The normalization is less 
homogeneous, but in directions well above the galactic disk and away from the galactic- 
centre region it has been found to be consistent with a normal distribution around the 
mean value: thus the claim of a possible extragalactic origin (Sreekumar et al. 1998). 

The origin of this diffuse GBR is still unknown. The published candidate sources range 
from the quite conventional to the decisively speculative. Perhaps the most conservative 
hypothesis is that the GBR is extragalactic, and originates from active galaxies (Bignami 
et al. 1979; Kazanas & Protheroe 1983; Stecker & Salamon 1996). The fact that blazars 
have a 7-ray spectrum with an average index 2.15 it 0.04, compatible with that of the 
GBR, supports this hypothesis (Chiang & Mukerjee 1998). The possibility has also been 
discussed that Geminga-type pulsars, expelled into the galactic halo by asymmetric super- 
nova explosions, be abundant enough to explain the GBR (Dixon et al. 1998; Hartmann 
1995). More exotic hypotheses include a baryon-symmetric universe (Stecker et al. 1971), 
now excluded (Cohen et al. 1998), primordial black hole evaporation (Page & Hawking 
1976; Hawking 1977), supermassive black holes formed at very high redshift (Gnedin & 
Ostriker 1992), annihilation of weakly interactive big-bang remnants (Silk & Srednicki 
1984; Rudaz & Stecker 1991), and a long etc. 

The EGRET GBR data in directions above the galactic disk and centre show a signif- 
icant deviation from isotropy, correlated with the structure of our galaxy and our position 
relative to its centre (Dar et al. 1999). In Fig. 4 we have plotted, as a function of 9, 
the EGRET GBR counting-rate above 100 MeV. This figure clearly shows, in three out 
of the four quarters of the celestial sphere, an increase of the counting rate towards the 
galactic centre. The "reduced" x^ per degree of freedom for a constant flux is 2.6: very 
unsatisfactory. A best fit of the form F = Fq + Fi (1 — cos 9) yields x^ = 1.3, a very large 
amelioration (for higher polynomials in cos 9 the higher-order coefficients are compatible 
with zero: the fit does not significantly improve). Note also that at angles with cos 9 larger 
than its mean value {cos 9) = 0.0246 (9 < 88.6°), 10 out of the 12 data points are above 



the average flux, while at angles with 9 > 88.6°, 18 out of the 24 data points are below 
the average. The probability for a uniform distribution to produce this large or a larger 
fluctuation is 1.5 x 10~^. Even in directions pointing to the galactic disk and the galactic 
centre, EGRET data on 7-rays above 1 GeV show an excess over the expectation from 
galactic cosmic-ray production of vr^'s (Pohl & Esposito 1998). Electron bremsstrahlung 
in gas is not the source of the 1-30 MeV inner-Galaxy 7-rays observed by COMPTEL 
(Strong et al. 1997), since their galactic latitude distribution is broader than that of the 
gas. These flndings also imply that inverse Compton scattering may be much more impor- 
tant than previously believed (Strong & Moskalenko 1998; Moskalenko and Strong, 2000; 
Dar et al. 1999). 

Thus, it appears that the EGRET data advocates a local (as opposed to cosmological) 
origin for the GBR. Indications of a large galactic contribution to the GBR at large 
latitudes were independently found by Dixon et al. (1998) by means of a wavelet-based 
"non-parametric" approach that makes no reference to a particular model. Strong & 
Moskalenko (1998) and Moskalenko & Strong (2000) also found that the contribution of 
inverse Compton scattering of galactic cosmic ray electrons to the diffuse 7-ray background 
is presumably much larger than previously thought. 

2. Galactic Origin of the GBR: In recent publications we went one step further 
(Dar et al. 1999; Dar &: De Rjula 2000) and explored in detail the possibility that the 
diffuse gamma-ray background radiation at high galactic latitudes could be dominated 
by inverse Compton scattering of cosmic ray (CR) electrons on the cosmic microwave 
background radiation and on starlight from our own galaxy. We made three very simple 
assumptions: (1) The locally-measured electron spectrum, shown in Fig. 4, which is well 
fitted, from E ~ 10 GeV to ~ 2 TeV by, 

E 



1-p r -c^ n-3.2±0.l0 
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cm"^ s"^ sr"^ MeV"\ (2) 



MeV 

is representative of its average form throughout the Galaxy and its extended CR halo. (2) 
Above a few GeV the acceleration of CR electrons is in equilibrium with their cooling by 
inverse Compton scattering (ICS) on starlight and on the microwave background radiation. 



(3) The GBR is produced by ICS of CR electrons on CBR and starlight photons in our 
Galaxy and in external galaxies. 

Assumptions (1) and (2) allow us to derive directly, the GBR spectral index from the 
observed electron index: (3^ = (/3c + l)/2 = 2.10 ± 0.05. It agrees with the GBR index as 
observed by EGRET, /3^ = 2.10 ± 0.03. This index is universal, independent of whether 
the ICS took place in our Galaxy or in external galaxies. 

The two dominant contributions to the GBR within our model are inverse Compton 
scattering of Galactic CR-electrons off the cosmic background radiation and starlight. 
There is a small additive effect of ICS from sunlight. The contribution from external 
galaxies, is also sub-dominant. The ICS spectrum -a cumbersome convolution (Felten 
&; Morrison 1966) of a CR power spectrum duc/dE = AE^'^^ with a photon thermal 
distribution- can be approximated very simply (Dar &: De Riijula 2000) by 
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where o"t = 0.65 x 10^^^ cm^^ is the Thomson cross section, ei are the mean photon 
energies of starlight and of the CBR and Ni(b, 1) are the effective CR electron column 
density resulting from the convolution of the space distribution of CR electrons with those 
of starlight and of the CBR in the direction (b,l). Simple analytical expressions for them 
are given in Dar & De Rujula (2000). 

Our predictions of the intensity and the angular dependence of the GBR are shown in 
Figs. 5, 6. Like that for the spectral index, they are in excellent agreement with those 
observed for the alleged extragalactic GBR, leading to the conclusion that the GBR could 
be dominated by the emission from our own galaxy. 

The only non-conventional aspect of our model is that, in order to reproduce the 
observed intensity of the GBR, we must assume the scale height of our galaxy's CR- 
electron distribution to be almost twice the traditionally-accepted upper limit. Such a 
large scale height is not in contradiction with radio synchrotron-emission from our galaxy 



if the galactic disk and its magnetic field are embedded in a larger magnetic halo with a 
much weaker field. 

3. Predictions: The main predictions specific to our scenario are: 

(1) The GBR should reflect the asymmetry of our off-centre position in the Galaxy. 

(2) Very nearby star-burst Galaxies, such as M82, and radio galaxies with presumably 
large CR production rates, such as Cygnus A, may be visible in gamma rays. 

(3) The GBR spectrum should not have the sharp cutoff, above E ~ 100 GeV, expected 
(Salamon & Stecker 1998) for cosmological sources. But it should nonetheless steepen 
around 10-100 GeV, because of the anticipated "knee" in the electron spectrum and of 
the energy-dependence of the Klein-Nishina cross section. 

These features of our scenario should be testable when the next generation of cosmic-ray 
and 7-ray satellites (AMS-02 and GLAST) are operational, hopefully by 2005. 
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Figure 1: Comparison between the spectrum of the GBR, measured by EGRET (Sreeku- 
mar et al. 1998), and the prediction for ICS of starlight and the CMB by CR electrons. 
The slope is our central prediction, the normalization is the one obtained for he = 20 kpc, 
Pe = 35 kpc. 
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Figure 2: EGRET data on the GBR spectral index as a function of 9 the angle away 
from the direction of the galactic center. The line is the predicted spectral index. The 
various plots correspond to the individual half- hemispheres, (a) b > 0, 1 > 0. (b) b > 0, 
1 < 0. (c) b < 0, 1 > 0. (d) b < 0, 1 < 0. 
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Figure 3: EGRET data, organized as in Fig. 2, for the dependence on 9 of the GBR 
intensity above 100 MeV. 
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Figure 4: The primary cosmic-ray electron spectrum (Evenson & Meyers 1984; Golden et 
al. 1994; Ferrando et al. 1996) as measured by Prince 1979 [crosses]; Nishimura et al. 1980 
[squares]; Tang 1984 [circles]; Golden et al. 1984 [triangles]; Barwick et al. 1998 [stars]. 
The slope is the prediction, the magnitude is normalized to the data. 
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Figure 5: Contributions to the GBR flux above 100 MeV as functions of longitude 1, at 
fixed latitude b, from ICS of starlight (dotted), and CMB in our galaxy (dashed); from 
the total ICS from external galaxies (continuous), and from sunlight (dot-dashed). The 
vertical scale is 10'^ times the number of photons/ [cm^ ssr]. The results are for he = 20 
kpc, pe = 35 kpc. 
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Figure 6: The flux of GBR pliotons above 100 MeV: comparison between EGRET data 
and our model for he = 20 kpc, pc = 35 kpc, as functions of latitude at various fixed 
longitudes. The grey domain is EGRET's mask. 
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